Observation of the Dynamical Casimir Effect in a 
Superconducting Circuit 

CM. Wilson^, G. Johansson^, A. Pourkabirian^, J.R. 

Johansson^, T. Duty^, F. Nori^'"^ & P. Delsing^ 
^Depart, of Microtechnology and Nanoscience, 
Chalmers University of Technology, Goteborg, Sweden 
^Advanced Science Institute, RIKEN, 
Wako-shi, Saitama 351-0198, Japan 
^University of New South Wales, Sydney, NSW, 2052 Australia and 
^University of Michigan, Ann Arbor, MI 48109, USA 
(Dated: January 25, 2013) 



One of the most surprising predictions of modern quantum theory is that the vacuum of 
space is not empty. In fact, quantum theory predicts that it teems with virtual particles 
flitting in and out of existence. While initially a curiosity, it was quickly realized that these 
vacuum fluctuations had measurable consequences, for instance producing the Lamb shift 
of atomic spectra and modifying the magnetic moment for the electron[2j. This type of 
renormalization due to vacuum fluctuations is now central to our understanding of nature. 
However, these effects provide indirect evidence for the existence of vacuum fluctuations. 
From early on, it was discussed if it might instead be possible to more directly observe 
the virtual particles that compose the quantum vacuum. 40 years ago, Moore [5] suggested 
that a mirror undergoing relativistic motion could convert virtual photons into directly 
observable real photons. This effect was later named the dynamical Casimir effect (DCE). 
Using a superconducting circuit, we have observed the DCE for the first time. The circuit 
consists of a coplanar transmission line with an electrical length that can be changed at a 
few percent of the speed of light. The length is changed by modulating the inductance of 
a superconducting quantum interference device (SQUID) at high frequencies (~ 11 GHz). 
In addition to observing the creation of real photons, we observe two-mode squeezing of the 
emitted radiation, which is a signature of the quantum character of the generation process. 

That a mirror can be used to measure vacuum fluctuations was first predicted by Casimir 
in 1948[1]. Casimir predicted that two mirrors, i.e. perfectly conducting metal plates, held 
parallel to each other in vacuum will experience an attractive force. Essentially, the mirrors 
reduce the density of electromagnetic modes between them. The vacuum radiation pressure 
between the plates is then less than the pressure outside, generating the force. As this static 
Casimir effect can then be explained by a mismatch of vacuum modes in space, the dynamical 
Casimir effect can be seen as arising from a mismatch of vacuum modes in time. As the 
mirror moves, it changes the spatial mode structure of the vacuum. If the mirror's velocity, 
V, is slow compared to the speed of light, c, the electromagnetic (EM) field can adiabatically 
adapt to the changes and no excitation occurs. If instead v/cis not negligible, then the field 
cannot adjust smoothly and can be nonadiabatically excited out of the vacuum. 

The static Casimir effect can also be calculated in terms of the electrical response of the 
mirrors to the EM field|5j. A similar complementary explanation exists for the DCE|3]. 
Theoretically, the ideal mirror represents a boundary condition for the EM field, in partic- 
ular, that the electric field is zero at the surface. This boundary condition is enforced by 
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the flow of screening currents in tlie metal. A mirror moving in a finite EM field then losses 
energy as the screening currents will emit EM radiation, as in an antenna. Classically, we 
expect this radiation damping to be zero in a field-free region. In quantum theory, however, 
the screening currents must always act against the vacuum fluctuations. Therefore, even 
moving in the vacuum will cause a mirror to emit real photons in response to the vacuum 
fluctuations. 

If we consider the literal experiment of moving a physical mirror near the speed light, 
we quickly see that this experiment is not feasible. Braggio et al. considered [0] the case of 
moving a typical microwave mirror in an oscillating motion at a frequency of 2 GHz with a 
displacement of 1 nm. This produces a velocity ratio of only f/c ~ 10^^ with an expected 
photon production rate of approximately 1 per day. Nevertheless, it requires an input of 
mechanical power of 100 MW while, at the same time, the system would need to be cooled 
to ~ 20 mK to ensure that the EM field is in its vacuum state. These difficulties have lead 
to a number of alternative proposals [71-[T^ . for instance utilizing surface acoustic waves, 
nanomechanical resonators, or by modulating the electrical properties of cavities. 

Here we investigate one such proposal using a superconducting circuit [I2l [13] : an open 
transmission line terminated by a SQUID. A SQUID is comprised of two Josephson junctions 
connected in parallel to form a loop. At the frequencies studied here, the SQUID acts as a 
parametric inductor whose value can be tuned by applying a magnetic flux, $ext; through 
the SQUID loop. When placed at the end of a transmission line, this SQUID can then be 
used to change the line's boundary condition. In previous work, we showed that this tuning 
can be done on very short time scales [HI [16]. The changing inductance can be described 
as a change in the electrical length of the transmission line and, in fact, provides the same 
time- dependent boundary condition as the idealized moving mirror [T7[ [T8] . In the same 
way as for the mirror, the boundary condition is enforced by screening currents that flow 
through the SQUID. Unlike the mirror, the effective velocity of the boundary, defined as the 
rate of change of the electrical length, can be very large, approaching f/c ~ 0.05 for a 10% 
modulation of the SQUID inductance. The photon production rate is therefore predicted to 
be several orders of magnitude larger than in other systems. 

Quantum theory allows us to make more detailed predictions than just that photons 
will simply be produced. If the boundary is driven sinusoidally at an angular frequency 
Ud = 2nfd, then it is predicted that photons will be produced in pairs such that their 
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frequencies, ooi and uj2, sum to the drive frequency, i.e., we expect ood = iOi + 002- This 
pairwise production imphes that the EM field at these frequencies, symmetric around a;d/2, 
should be correlated. In detail, we can predict that the field should exhibit what is known 
as two-mode squeezing [19j. These correlations are a signature of the two-photon nature of 
the photon generation process. 

Theoretically, we treat the problem as a scattering problem in the context of quantum 
network theory |20]. For superconducting circuits, it is convenient to describe the EM field 
in the transmission line in terms of the phase field operator t) = E(x, t')dt', where 
E{x,t) is the electric field operator. In the CPW, is described by the massless, 

scalar Klein-Gordon equation in one dimension, the solution of which can be written as 
t) = 0in(x — cot) + 0out(a; + cot), whcrc 0in(out) is the field propagating inward to (outward 
from) the SQUID and cq ~ 0.4c is the speed of light in the transmission line. We solve the 
scattering problem in Fourier space defining 




where a{u) and its hermitian conjugate a'^lu) and are the standard annihilation and creation 
operators and k^^ = u/cq is the wavenumber of the radiation. Solving the scattering problem 
then amounts to finding expressions for al^^{u) and aout(w) as a function of al^^u) and ain(w). 
The boundary condition imposed by the SQUID determines the connection between these 
operators. With the output operators, we can then calculate the properties of the measurable 
output field assuming the input field is in a definite state, such as a thermal state or vacuum 
state. For a static magnetic flux, $cxt7 we obtain the simple expressions aout(i^) = -R(w)ain(w) 
where R{uj) is the reflection coefficient from the SQUID. R{uj) = — exp[2iku]ie{^ext)] has 
the simple form [12] of a phase shift due to a transmission line of fixed length iei^e^t) = 
Lj{^cxt)/ Lq. Here Cq = 1/\/LqCq, Lq (Cq) is the inductance (capacitance) per unit length 
of the line, Lj($ext) = ('^'o/27r)^/i?j(<l'cxt) is the Josephson inductance of the SQUID, Ej is 
its Josephson energy, and $0 = ^/2e is the superconducting flux quantum. 

In order to generate DCE radiation, must change with a nonuniform acceleration. A 
simple example of this type of motion is just a sinusoidal drive with an amplitude of Sie- If 
it is driven at Ud, we then find a simple expression [13] for aout(i^) in the region u < u^- 

a-ontiuj) = R{uj)ain{uj) + S{uj)alX^d - c^) (1) 
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where S{uj) = —i^Sle/coj^/^ij^d — oo)A{uj)A*{iOd — uj) and A{uj) is the spectral amphtude of 
the transmission hne. Crucially, the time-dependent boundary leads to mixing of the input 
field's creation and annihilation operators. With this expression we can calculate the output 
photon flux density for an input thermal state 

nout(w) = {al^^{u)aont{uj)) = riiniu) + \S{u)\^nin{ud -uj) + (2) 

The first two terms, proportional to ni^{uj), represent the purely classical effects of reflection 
and upconversion of the input field to the drive frequency. They are zero at zero temperature. 
The last term is due to vacuum fluctuations and is, in fact, the DCE radiation. 

The photon production rate depends on the density of photonic states in the transmission 
line, which is |yl(a;)p. For an ideal transmission line, A{u) = 1 and the DCE radiation 
is n^^^{ijj) = {Sie/coY'^i^d — <-^)- The integrated photon flux of the DCE radiation is 
Edge = (ij-'d/127r)(t'e/co)^ where Ve = Sie'^d is the maximum velocity of the boundary. The 
relativistic nature of the effect is apparent here in that the photon flux goes to zero if we 
allow the speed of light to go to infinity. Finally, we note that this spectrum is identical to 
that calculated for an ideal mirror oscillating in 1-D space [T8]. 

We present measurements on a sample with a short (~ 100 fim) Al coplanar waveguide 
(CPW) on-chip which transitions to a Cu CPW on a microwave circuit board (see Fig. la 
& b). As a first measurement (see Fig. Id), we can measure the phase shift of a microwave 
probe signal reflected from the SQUID as we change $ext- This illustrates how the SQUID 
changes the boundary condition at the end of the line. 

We can now look at the effects of nonadiabatic perturbations. The flux through the 
SQUID is driven at microwave frequencies by an inductively coupled CPW line that is 
short-circuited ~ 20 fim from the SQUID. The sample is cooled to 50 mK in a dilution 
refrigerator. This corresponds to a thermal photon occupation number of n = 0.008 at 5 
GHz, the center of our analysis band. If we consider the last two terms in Eq ^ , which 
are the response of the system to the changing boundary, we can compare this small value 
n = 0.008 to the vacuum response, which has a coefficient of 1. We therefore conclude 
that thermal effects are negligible and disregard them in the rest of the paper. For the 
first measurements, the drive is operated in continuous wave (CW) mode and the analysis 
frequency of the digitizer tracks the drive at Ud/2 with a 100 kHz bandwidth determined by 
digital filtering. We study the total output power as a function of drive frequency and power. 
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FIG. 1: a) Optical micrograph of the device. Light parts are Al while dark parts are the Si 
substrate. The output line is labeled "CPW" and the drive line enters from the top. Both lines 
converge near the SQUID, b) A scanning-electron micrograph of the SQUID. The SQUID has a 
normal state resistance of 218 0, implying a Josephson inductance at zero field of Lj{0) = 0.23 
nH. c) A simplified schematic of the measurement setup. In addition to the driving line, a small 
external coil is used to apply a dc flux bias. The driving line has 36 dB of cold attenuation along 
with an 8.4 to 12 GHz bandpass filter. The filter ensures that no thermal radiation couples to the 
transmission line in the frequency region were we expect DCE radiation. The outgoing field of 
the CPW is coupled through two circulators to a cryogenic HEMT amplifier (LNA) with a system 
noise temperature of T/v ~ 6 K. The signal is further amplified at room temperature before being 
captured by two vector RF digitizers which use heterodyne mixing followed by digitization of the 
intermediate frequency signal, d) Measured phase shift of a microwave probe signal reflected from 
the SQUID, as it is tuned with a static flux. The probe frequency is 6.18 GHz. 
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The results are shown in Fig. 2a. We clearly see photon generation for essentially all drive 
frequencies spanning the 8-12 GHz band set by the filtering of the line. This corresponds to 
an analysis band of 4-6 GHz. The detailed drive power dependence of the output depends 
on the frequency and reflects a combination of the frequency response of both the drive line 
and the output line. 

We can quantify the photon production rate by comparing it to the amplifier noise tem- 
perature. We see that the produced photons roughly double the noise level, suggesting a 
power per unit bandwidth of a few Kelvin. We predict that in an ideal transmission line that 
the power should instead be a few mK. [T2] In general though, transmission lines are not ideal 
and have parasitic resonances with a low quality factor {Q) associated with reflections at 
connectors, etc.. These resonances modify the electromagnetic density of states in the trans- 
mission line, |74(ci;)p, thereby enhancing[T3] Fdce by a factor of ~ Q^. In our measurement 
setup, we see such resonances, with Q ~ 30-50, likely associated with the superconducting 
cable which runs between the sample and the amplifier. This implies an enhancement of the 
photon production rate by a factor of 1000-2000, consistent with what we observe. 

In the next set of measurements, we fix the drive frequency, but scan the digitizer fre- 
quency. In this way, we can see over what band photons are produced for fixed drive 
frequency. In Fig. 2 b & c, we show the results for two different drive frequencies. We 
clearly see broadband photon production at all analysis frequencies, including detunings 
from ffi/2 larger than 2 GHz. 

The broadband nature of the photon generation, both in terms of the drive frequency and 
analysis frequency, clearly distinguishes the observed phenomenon from that of a parametric 
amplifier comprised of a driven single-mode oscillator or cavity. For a parametric amplifier, 
we expect the photon production to be narrow band in both drive and analysis frequencies, 
which is clearly not the case in this experiment. 

Theory also predicts that the output should exhibit voltage-voltage correlations at differ- 
ent frequencies with a particular structure commonly known as two-mode squeezing (TMS). 
Following Ref. |19]. we can describe a two-photon state, as we expect the DCE to gen- 
erate, in terms of modulation of the center frequency of the state. We can then define 
the modulation operators at the detuning e as ai(e) = (A+a(w-|.) -|- X^a^{ijj^ ))/V2 and 
Q^2(e) = {—i\+a(uj^) + iX^a){ujJ))/ \/2 where IjJ± = OJ^jl ± e and \± = {2ijj±/ijJdY^'^ ■ The 
factors \± rescale the operators from quanta at uj± to quanta at the center frequency u;d/2. 
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FIG. 2: Photons generated by the dynamical Casimir effect. Here we show the output power of the 
transmission hne while driving the SQUID. Black is low power, while yellow is high power, a) The 
drive frequency and power are scanned in continuous wave mode . The digitizer tracks the drive 
frequency at fd/2 with a 100 kHz analysis bandwidth, b & c) Broadband photon generation. The 
drive frequency is fixed while the digitizer frequency is scanned. The drive is chopped on and off 
and we record the difference in the power level for the on and off state. The on (off) time is 50 ms 
(50 ms) and we measure in a 50 MHz analysis bandwidth. The drive frequencies in b) and c) are 
8.70 GHz and 11.30 GHz respectively. Half the drive frequency is indicated by the blue line. We 
see that for the two drive frequencies, the spectral shape of the output is similar. For instance, the 
output is more intense between 4-4.5 GHz in both cases. This shows how the spectral density of 
the output line affects the radiation. (The decline of the power below 4 GHz is due to the roll-off 
of the amplifier gain.) We also see that when half the drive frequency corresponds to the output 
mode, the necessary drive power decreases. 

We see that these operators mix excitations at the upper and lower sidebands of the field 
with a definite phase. The TMS of the field then appears as an imbalance of the noise 
in one of these modes compared to the other. We can then define a normalized squeezing 
statistic pi] 

Ell — 



^11 — ^22 /„N 

^2 = y , y (3) 
2^11 + 2^22 

where S^n = {ctmO:}^ + al^am)/2 is the symmetrized spectral density matrix. For the case of 
an ideal transmission line, we can calculate the TMS of the output field to be 

,uj+uj- 6ip 1 



^2 = 2-^— — = . (4) 

Ud Co 1 + w+w_(54/co)^ 

For small detunings (e ^ uJd/2), we get the approximate expression a2 ~ fe/2co, which 
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implies a squeezing of about 2% for a 10% modulation of the SQUID inductance. 

Experimentally, we measure the four quadrature voltages of the upper and lower side- 
bands I± and Q±. The observable (hermitian) quadrature operators can be related to 
creation and annihilation operators as 

I± = {a{^±) + a{^±)^) ; Q± = -^\l^^^^{a{^±) - a{^±)^)- (5) 

We can then write o"2 in terms of the quadratures as 

^2 = ((/+/-) -(g+g-)). (6) 

where Pavg = + {Q+) + (-^^ ) + {Q-))/'^ is the average noise power in the sidebands. The 
theory also predicts a special structure for the correlations of the TMS, in particular that 
(/+/_) = —{Q+Q-) and that = {I-Q^). Finally, we comment that S transforms 

under phase rotations (see eqn. ([T])) such that we can specify (/+g_) = (/_g+) = without 
loss of generality, which has been done in writing eq. ([6]). 

To measure the correlations, we use a single amplifier but take advantage of the fact 
that amplifier noise at different frequencies is uncorrelated. After amplifying, we split the 
signal and feed the two outputs into two separate digitizers which are synchronized. We 
then calculate the four IQ cross-correlation functions. [22j Typical results are shown in Fig. 
3. In Fig. 3a & b, we show the (/+/-) and {Q+Q-) cross correlations with the drive on 
and drive off for comparison. With the drive on, we see very clear cross correlations which 
are ~ 100 times larger than the parasitic amplifier correlation. Also, we see that indeed 
(/+/_) = — (g+g_) as we expect for two-mode squeezing. We also see that this is not the 
case for the amplifier noise. The correlations imply a value of cr2 ~ 0.04 consistent with 
our expectations. The data in Fig. 3 is measured with a frequency seperation of 40 MHz 
(e/27r = 20 MHz), although we have measured similar values with separations as large as 
700 MHz. 

Theory predicts that even as the field is two-mode squeezed, if we look at either sideband 
frequency individually, we expect it to remain unsqueezed, essentially appearing as a thermal 
field at some effective temperature. In Fig. 3c, we then plot the TMS of the field, (T2, along 
with the one-mode squeezing, ai = ((/^) — {Q^))/{{P) + {Q^)) at both separate frequencies, 

and W-, as a function of drive power. Although there is some scatter in (T2, we see that it 
clearly increases as a function of drive power while the one-mode fields remain unsqueezed. 
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FIG. 3: Two-mode squeezing of the DCE field, a) and b) The normalized cross-correlation functions 
(/+/_) /Pavg and {Q+Q-)/Pe,vg with the drive on (a) and drive ofi' (b). Note that the scale of (b) 
is expanded by xlOO. The drive frequency is 11.30 GHz and the two digitizer frequencies are 
5.63 GHz and 5.67 GHz. The RF bandwidth is 10 MHz and 4 x 10'' samples of each voltage 
are used. With the drive on, (a), we clearly see cross correlations of a few percent. Note that 
(/_!-/_) = —{Q+Q-) as predicted. (The shape of the correlation functions in delay time is largely 
determined by the antialiasing filters of the digitizers.) With the drive off, (b), we see that the 
residual correlations of the amplifier are negligible, c) The two- mode squeezing, cr2, of the field 
along with the one- mode squeezing, cji, at both a;+ and uj- as a function of drive power. Although 
there is some scatter in cj2, we see that it clearly increases while the single-mode fields remain 
unsqueezed, d) The effects of phase rotations on the complex correlation function ^ (defined in 
the text). The color scale is the real part of ^. The x-axis corresponds to a change in the drive 
phase. The y-axis corresponds to a digital rotation of the measured value of ^. As expected for 
TMS, we see that ^ is symmetric under phase rotations. 
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This is an important additional check that the correlations arise from two-mode squeezing, 
and not just some well-timed modulation of, for instance, the amplifier gain that would 
produce the same type of apparent squeezing at the offset frequencies. 

For a two-mode squeezed state, we can predict how the correlations transform under 
rotations of the phase of the EM field by an angle 9. In particular, if we define the appropriate 
combination of correlation functions 

^ = ((/+/_)) - (g+g_)) + t{{i+Q-) + (/-g+)) (?) 

we expect ^ to transform such that ^' — e~^*^^. To explore this predicted symmetry, we 
can compute the complex quantity ^ from the experimental correlation functions and look 
at the rotation properties. In Fig. 3d, wc compare the results of rotating the measured 
quantity \1' to changing the drive phase. We see, as expected for a TMS state, that the 
rotation of one cancels a rotation of the other. 

In conclusion, we observe broadband generation of microwave photons in an open trans- 
mission line with a periodically modulated boundary condition. The emitted photons exhibit 
two-mode squeezing correlations, which are characteristic of photons generated in correlated 
pairs. Taken together, wc believe these results represent the first experimental observation 
of the dynamical Casimir effect. 
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